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: ABSTRACT 

Collective flavor oscillations driven by neutrino- neutrino interactions inside core-collapse super- 
novae have now been shown to drastically alter the resultant neutrino fluxes. This would in turn 
. signiflcantly affect the diffuse supernova neutrino background (DSNB), created by all core-collapse 
' supernovae that have exploded in the past. In view of these collective effects, we re- analyze the 
[ potential of detecting the DSNB in currently running and planned large-scale detectors meant for 
■ detecting both i/g and Ug. We flnd that the event rate can be different from previous estimates by 
' upto 50%, depending on the value of 6*13. The next generation detectors should be able to observe 
DSNB fluxes. Under certain conducive conditions, one could learn about neutrino parameters. 
For instance, it might be possible to determine the neutrino mass hierarchy, even if 6*13 —>■ 0. 
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1 Introduction 



Observation of neutrino signal from a core-collapse supernova (SN) is expected to contribute 
significantly towards determination of the neutrino mass hierarchy and the mixing angle ^13. 
It is also expected to shed light on SN shockwave dynamics and set stringent bounds on the 
existence of sterile neutrinos. Thus a lot of attention has been devoted to this scenario of resonant 
neutrino conversions in a SN [1, 2]. The work has focussed on the determination of mass hierarchy 
and signatures of a non-zero ^13 [3, 4], Earth matter effects on the neutrino fluxes when they 
pass through matter [5, 6, 7], shock wave effects on observable neutrino spectra and their model 
independent signatures [8]- [13]. Recently, possible interference effects for multiple resonances [14], 
sterile neutrinos [15], the role of turbulence in washing out shock wave effects [16, 17, 18], and time 
variation of the signal [19] have also been explored. Observational data on SN neutrinos come 
only from the SN1987A [20, 21, 22], and while it confirms broad expectations, the small number 
of events limits our ability to draw detailed inferences. 

However, supernovae are relatively rare in our galaxy with an estimated rate of about 1 — 3 
per century [23], which prompts consideration of the alternative strategy to detect neutrinos 
from supernovae that are further away. Neutrinos accumulated in the Universe from all the 
SN explosions in the past and present epoch form a cosmic background, known as the diffuse 
supernova neutrino background (DSNB) or supernova relic neutrinos [24, 25, 26]. The expected 
flux of these DSNB neutrinos depends mainly on the SN rate and the "typical" flavor dependent 
flux of neutrinos from supernovae. 

The SN rate can be either determined directly [27] or from the cosmic star formation rate, 
which is measured using a variety of ways like galaxy luminosity function of rest-frame ultravi- 
olet radiation [28]-[36], far-infrared/sub- millimeter dust emission [37, 38] and near-infrared Ha 
fluorescent line emission [39]- [42] and radio emission [43]. Though these techniques suffer from 
various ambiguities and complications like dust extinction [44], careful studies have enabled a 
precise determination of the star formation rate [45] . 

The typical neutrino flux from a SN on the other hand is not experimentally available. The 
data from SN1987A do not allow a clean determination of the spectral parameters [46]. So one 
has to resort to using primary neutrino fluxes predicted using SN simulations. Fluxes predicted 
by different groups (and sometimes different simulations by the same group) are at considerable 
variance, because of their different physics input [47, 48]. Additionally, these primary fluxes are 
further mixed by flavor conversions as they stream through the SN, thus requiring knowledge 
of not only the typical primary spectra, but also the typical SN density profiles. The emitted 
neutrino fluxes are therefore ridden with uncertainities at present. However, they could be made 
precise with more sophisticated simulations or observation of a galactic SN. 

Estimation of the DSNB flux has been performed in previous literature, with varying ap- 
proaches and results [45], [49]-[62]. Most studies have focussed on DSNB detection via Ve scat- 
tering off protons at water Cerenkov detectors [61] and large liquid scintillator detectors [55]. On 
the other hand, Vf, detection has been considered at a liquid argon detector [63] and at Sudbury 
Neutrino Observatory (SNO) [64, 65]. In [62], authors have performed a detailed comparative 
study of Vf. detection in different future large scale observatories - by interaction of on oxy- 
gen in water Cerenkov detectors, on carbon in liquid scintillator detectors and on argon in liquid 
argon detectors. Experimentally, the best upper hmits at 90 % C.L. of 6.8 x 10^ v^. cm~^s~^ 
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(25 MeV< E^^ < 50 MeV) and 1.2 iy^ cm-^s^i {Ep^ > 19.3 MeV) come from the Liquid Scintil- 
lation Detector (LSD) [66] and the Super-Kamiokande (SK) detectors [67] respectively. However, 
stronger bounds can be placed on these fluxes, albeit using somewhat indirect arguments [65, 68]. 
Some of the theoretical estimates of the DSNB fluxes predict event-rates for z/g that are tanta- 
lizingly close to detection, e.g., the observational upper limit set by the SK collaboration [67]. 
The prospects for discovery thus seem promising if a large water Cerenkov detector like SK is 
loaded with 0.02% GdCls [69] or if one or more of the proposed next generation detectors become 
available. 

The study presented in this work removes a common assumption, made for all previous esti- 
mates, that neutrino-neutrino interactions are too feeble to cause any flavor conversion. Although 
these interactions had been studied in previous literature [70]-[75], it has only very recently been 
appreciated that they induce sizable flavor conversion in supernovae [76, 77]. Motivated by this 
interesting result, the effects of these neutrino-neutrino interactions have been explored in the 
context of supernovae, in a scries of papers [78]- [92]. The interesting aspect of these conversions 
is that the neutrinos and antincutrinos of different energies imdergo conversions together, and are 
almost in-phase. Therefore these conversions are referred to as being "collective". The effect of 
these interactions becomes negligible beyond the first few hundred kilometers, when the neutrino 
densities become much lower. The end of these collective effects is marked by a complete (or 
step-wise) swapping between the flavor spectra of antincutrinos (or neutrinos) for inverted hier- 
archy (IH) [80, 82, 83]. Further out from the centre of the star, the traditional picture of flavor 
evolution by Mikheyev-Smirnov-Wolfenstein (MSW) conversion is not changed, except that the 
primary fluxes emitted at the neutrinosphere undergo the above-mentioned "pre-processing" due 
to the collective effects. The fluxes emitted by a SN have already been calculated in a three-flavor 
framework, including the effect of collective oscillations [89]. In this paper, we take these SN neu- 
trino fluxes calculated in [89] , the SN rate deduced from the cosmic star formation rate calculated 
by Beacom et al. [45], and the standard A-CDM cosmological model [93] as inputs to calculate 
the DSNB flux. The expected DSNB flux in the case of IH turns out to be quite different from 
those contained in previous works that disregarded collective effects. Thus the prospects of DSNB 
detection at antineutrino and/or neutrino detectors are modified. We report the DSNB fluxes 
and their observability, with and without neutron tagging, at present and proposed detectors. 
These changed expectations have impact on limits that can be set on non-standard models and 
interactions of neutrinos as well. 

The paper is organised as follows. In Section 2, we outhne the dependence on the DSNB event- 
rate at Earth as a function of the cosmic star formation rate and the primary neutrino fluxes and 
mention our choice of models for the same. We discuss our choice of detectors. In Section 3, we 
calculate the event-rates, dependent on the above choices, and present the results. We conclude 
the paper in Section 4, by summarizing our results and giving an outlook on the impact of these 
results. 

2 DSNB in Terrestrial Detectors 

From the early times to the present date, SN explosions have been fairly common events in the 
Universe. These explosions have injected a large number of neutrinos, with energies of tens of 
MeV, in the Universe. These neutrinos have created a diffuse background of SN neutrinos known 
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as diffuse supernova neutrino background. Evidently, tlie DSNB ffux depends on two ingredients: 

• The rate of SN explosions Rsn{z), as a function of cosmological redshift z. 

• The differential flux of neutrinos F^{Ei,), from a typical core-collapse event at redshift z. 

The differential flux of neutrinos Fi,{E^) depends on the primary neutrino fluxes F^{Ei,), emitted 
from the neutrinosphere, which get modified due to 

• Collective effects, i.e. neutrino-neutrino self interaction, close to the neutrinosphere. 

• MSW effects, i.e matter driven neutrino oscillations in the SN mantle and envelope. 

The total differential DSNB flux arriving at terrestrial detectors, expressed as the number of 
neutrinos of flavor u (where v — i^e, Vt and antineutrinos are denoted with a bar overhead) 
arriving per unit area per unit time per unit energy, due to all supernovae in the Universe up to 
a maximum redshift Zmaxi is 

/■o dt 

Fl{E,) = / {dz —) (1 + z) Rsn{z) + z)E,) . (1) 

Zmax ^Z 

Here E^, is the neutrino energy at Earth and Rsn{z) is the SN rate per comoving volume at 
redshift z. For our numerical calculations we have assumed Zmax — 7. Note that the factor (l + z) 
in the neutrino spectrum + z)Eiy) incorporates the redshift of the energy spectrum. 

From the Friedmann equation for a flat universe we have 

dy 

^ = -Hoil + + zf + n^y/' . (2) 

Thus the differential number flux of DSNB is 

Ho Jo + zy + Qa) 

For the standard A-CDM cosmology, we have 

fi„ = 0.3 ; = 0.7 and Hq = 70 hjo km s"^ Mpc"^ . (4) 

Therefore, we only need to know the SN rate Rsn{z) and the differential flux of neutrinos Fi,[Ei,), 
from a typical core-collapse event to calculate the DSNB flux at Earth. 

2.1 The Cosmic Supernova Rate 

The SN rate Rsn{z) is related to Rsf{z), through the initial mass function f{m), which describes 
the differential mass distribution of stars at formation [45, 56]. We assume that all stars that are 
more massive than SM© give rise to core-collapse events and die on a timescale much shorter than 
the Hubble time, and that the initial mass function (f{m) is independent of redshift. This allows 
us to relate the star formation rate Rsf{z) to the cosmic SN rate Rsn{z) as 

Rsn{z) = Rsf{z) ° . (5) 

Io.iMq ^{m)mdm 
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For our estimates, we use the initial mass function from reference [96], i.e. 



, , f m-i-^o (O.IM0 < m < O.5M0) 

^ [ m ^-^^ (m > O.5M0) ^ ^ 

Putting the above expression into Eq. (5) we find 

Rsn{z) = 0.0132 Rsf{z)M-' . (7) 

It should be noted that the factor connecting Rsn and Rsf is quite insensitive to the upper limit 
of the integrations in Eq. (5). 

Recent careful studies on different indicators of the cosmic star formation rate have been used 
to calculate the Rsf and its normalization. We use the cosmic star formation rate per comoving 
volume, Rsf, from the concordance model advocated in [94, 95], which is given by 



. (1 + z < 0.97 

Rsf{z)(x\ + 0.97<z<4.48 , (8) 

^ (1 + ^)-^-^ 4.48 < z 



with the local star formation rate given by 

Rsf{Q) = 0.0197 Moyr-^Mpc"^ . (9) 

This model satisfies the experimental upper limit on DSNB set by SK [67], and hence is known as 
the concordance model [45, 97]. 



2.2 Neutrino Fluxes from Core-collapse Supernovae 
2.2.1 Primciry Neutrino Fluxes 

A core-collapse SN is typified by conversion of its gravitational binding energy of about 10^^ ergs 
to neutrinos and antineutrinos with energies of tens of MeV. This translates to a total flux of 
more than 10^^ per SN explosion. The initial luminosity is about 10^^ ergs/sec (emitted purely 
as Ue in the neutronization burst lasting for about 25 msecs), and slowly reduces to about 10^^ 
ergs/sec over the following ~ 10 sees. Since the event-rate for DSNB is not very sensitive to the 
neutronization burst, we can ignore it for our analysis ^. Subsequently, after the neutronization 
burst, neutrinos and antineutrinos of all three flavors are emitted with a pinched thermal spectrum, 
that is conveniently parametrized as [98] 

where L° is the luminosity in the flavor u, {E^) is the average energy of i/, and Ci/ is the pinching 
parameter at the neutrinosphere. As a notational convenience, since there is no difference expected 
between the // and r flavors for SN neutrinos, we will refer to them together as x. 

^The neutrinos in the burst phase undergo eompHcated flavor eonversions, particularly for a class of low-density 
SN with degenerate 0-IMg-Nc> eores [86, 88, 90, 91]. However, these issues are not likely to be important for DSNB 
due to the relative smallness of the integrated flux from the burst-phase. 
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Figure 1: The DSNB flux spectrum arriving at Eartli as a function of tlie (anti)neutrino energy at 
Eartli. Tlie upper panels show the Ve flux while the lower panels show the v^. flux. The left panels 
correspond to the hypothetical case where we have only MSW matter effects in the SN while the 
right panels correspond to the case where we have both collective as well as MSW-driven flavor 
transitions. We have assumed that the initial neutrinos are given by the Gl model. 
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For our study, we take 3 sets of representative values of $° = VI/ {E^), {Ey) and motivated 
by SN simulations. One simulation by the Lawrence Livermore group (LL), and two different 
simulations by the Garching group (Gl, G2) have been chosen for our estimates, as shown in 
Table 1. 



Model 








$0 




LL 


12 


15 


24 


2.0 


1.6 


Gl 


12 


15 


18 


0.8 


0.8 


G2 


12 


15 


15 


0.5 


0.5 



Table 1: The parameters of the used primary neutrino spectra models motivated from SN simula- 
tions of the Garching (Gl, G2) and the Lawrence Livermore (LL) group. We assume (p^ — 4 and 
= 3. The total luminosity is chosen to be 3 x 10^^ erg. 

Note that the LL simulation obtained a large hierarchy {E^^) < (Ep^) < {E^^) ^ {EpJ, and 
an almost complete equipartition of energies between flavors. The Garching simulations predict a 
smaller hierarchy between the average energies, incomplete equipartition, and increased spectral 
pinching. The differences in the values of these parameters arise from the different physics inputs. 

2.2.2 Collective Effects and MSW Transitions 

The primary fluxes (at the neutrinosphere) are further processed by collective effects and MSW 
conversions before they get emitted from the SN^. 

Near the neutrinosphere, due to the large neutrino density, the neutrino-neutrino interaction 
energy is very large. This ensures that the neutrinos exhibit synchronized oscillations, i.e. neutri- 
nos of all energies oscillate coherently with the average frequency. These oscillations do not give 
rise to any effective flavor conversion since the effective mixing angle is highly supressed due to the 
large MSW potential. As the neutrinos stream outward, the neutrino density becomes smaller, 
and bipolar oscillations begin to take place. In the case of IH, these oscillations have large am- 
plitude even for a vanishingly small mixing angle. These oscillations thus can lead to a complete 
swapping of the Pg ^-nd spectra. The and spectra cannot swap completely, because of 
lepton number conservation, and the swap occurs only above a certain energy Ec, giving rise to a 
spectral split [80]. Eventually, beyond a few hundred kilometers, the neutrino- neutrino interaction 
energy becomes negligible, and collective effects cease to be important. 

Thus for normal hierarchy (NH), the collective effects do not affect the fluxes significantly and 
only MSW conversions are at work. In particular, the MSW resonances affect the z/g flux, while 
the i?e flux remains almost unaffected. For IH, the collective effects swap the Ue and above a 
certain energy Ec, determined by lepton number conservation [80, 90, 91]. Assuming = F^^, 

^The detailed picture of collective effects presented herein is valid only for initial spectra that resemble the LL 
model. However we are interested in seeing the maximum effect that these new effects can cause, and for that 
purpose it suffices to ignore more complicated features in the spectrum [99] for Gl or G2 like models. 
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the split-energy Ec is given in the adiabatic approximation by the imphcit equation 

Jj' dE (f° (E) - = f dE {Fl{E) - Fl{E)) . (11) 

On the other hand for antineutrinos, all z/g and are swapped. This pre-processed flux now 
undergoes the traditional MSW conversions which now affect the z/g flux, and not the z/g flux. 
The neutrinos then travel independently (while getting redshifted) as mass-eigenstates until they 
reach Earth, wherein they are detected as flavor eigenstates before or after having undergone 
regeneration inside the Earth. The fluxes of ^-nd Pg arriving at Earth are given in Table 2, as 
calculated in [89]. The quantities such as , are the initial SN neutrino fluxes while F^^ are 
the resultant fluxes emerging from the SN at redshift z. The quantities stand for 



Normal hierarchy 




Inverted hierarchy 




F.. = si, [P,,Fl + (1 - P,,)Fl) 


+ cl,Fl 


p _ \ slA + cl,Fl 


{E < E,) 
{E > E,) 


F.. = cl,Fl + sl,Fl 




^Pe = Sjjl + Cl, ((1 - 


Pi.)Fl + ) 



Table 2: Electron neutrino and antineutrino spectra emerging from a SN. See [89] for a prescription 
for calculating these final spectra in terms of the primary spectra. 



sin^ 612 (taken to be 0.3 for numerical studies) and cos^ 9i2 respectively and P13 is the effective 
jump probability between the neutrino mass eigenstates due to the MSW resonance(s), and takes 
a value between and 1. The value of P13 is approximately if ^13 is large (i.e. 61^ ^ 6 degrees) 
and for smaller values of ^13 it has a non-trivial dependence on energy and time, due to multiple 
resonances [14, 19] and turbulence [16, 17, 18]. However, due to the small number of events, 
we can probably neglect these sub-leading effects that occur in the small time-window when the 
Shockwave is in the resonance region. 

To calculate the DSNB flux at Earth F'{E,y), we need to integrate the fluxes in Table 2, correctly 
redshifted and weighted by the SN rate Rsn{z), over redshift z. We show in Fig. 1 the DSNB Pe 
(upper panels) and i^e (lower panels) fluxes arriving on Earth as a function of their (anti)neutrino 
energy at Earth. Wc have assumed the Gl model for generating this flgure. Note that the energy 
spectrum gets degarded to smaller energies due to redshift. The left panels show the predicted 
fluxes when one takes both collective as well as MSW transitions into account. To bring out the 
contrast with what the situation was earlier, we show in the right panels the predicted fluxes if 
one does not take collective effects. Wc can see that for NH the prediction have remained the 
same even after collective effects were taken, whereas for IH the fluxes are completely different. 

2.3 Terrestrial Detectors 

An array of existing and planned detectors could catch the DSNB neutrinos. In what follows, we 
will consider in particular three types of detectors for observing DSNB Pe' 

• Water Cerenkov detectors 

• Liquid scintillator detectors 
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• Gadolinium loaded water Cerenkov detectors 

Detection of z/g is more difficult. Both water and liquid scintillator detectors can in principle detect 
(as well muon and tau flavored neutrinos and antineutrinos) . In water Cerenkov detectors this 
can be done through neutrino-electron scattering. On the other hand, in liquid scintillators in 
addition to the neutrino-electron scattering, one can detect through charged current interac- 
tion on ^^C, while the other species can be detected through the neutral current interaction on 
^"^C. However, the cross-section for these processes are rather low. Another detector technology 
that has been proposed for detecting i/g is to use a high Z material, such as lead (and/or iron), 
interleaved with scintillators. Among such proposals are the OMNIS/ ADONIS projects and the 
HALO experiment at SNOLAB [100]. Therefore the only chance for detecting the Ve DSNB would 
be in a reasonably large 

• Liquid argon detector 

2.3.1 Water Cerenkov Detectors 

An upper bound on the DSNB flux already exists from non-observation of these neutrinos at the 

SK experiment [67]. Using 1496 days of data with 22.5 kton of fiducial volume, the DSNB fiux 
has been constrained to be less than 1.2 cm^^ s^^ for 19.3 MeV < -Ej^ < 30 MeV. SK is still 
running and could provide further constraint or evidence for DSNB fluxes in the future. Megaton 
water detectors with flducial volume in the ballpark of 500 kton have been planned in Japan 
(Hyper-Kamiokande (HK)) [101], Europe (MEMPHYS) [102], and USA (UNO) [103]. These have 
been proposed to serve as the far detector for long baseline experiments with powerful accelerator 
beams. At the same time, they would be used to study neutrinos from natural sources, such as the 
Sun, atmosphere and nearby supernovae. In particular, they will be useful tools for the observation 
of DSNB fluxes. While in principle water detectors can detect neutrinos and antineutrinos of all 
flavors, the easiest to observe is Pg, which is captured on protons via the inverse beta decay process 

Ve+P^e'^ + n. (12) 

The emitted positron is observed through the Cerenkov cone produced by it. The "true" positron 
energy is approximately related to the neutrino energy by E^, — 1.3 MeV. The other types of 
neutrino species would scatter electrons and thereby could also be detected. However, the cross- 
section for neutrino-electron scattering is much lower compared to the reaction (12). Therefore, 
in this paper we will consider the detection of only Pg in water Cerenkov detectors. The number 
of events per kton of detector mass is given as 



utT dEj dE, Fl{E,)(j{E,)R{E,,E^) , (13) 



e 



where ut is the number of protons in a kton of detector mass, T is the total exposure time. Eg the 
measured positron energy, i?'""" is the lower energy threshold, E'^^ is the upper energy threshold, 
Fl,{E^) is the DSNB flux at Earth, a{E^) is the cross-section and R{Ei, — 1.3, E^) is the energy 
resolution of the detector. For the energy resolution we assume a Gaussian form 



^2 
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where all quantities are given in units of MeV and is the half width at half maximum (HWHM). 
For the water Cerenkov detector we use 

ge = 0.47^(S^(MeV) - 1.3) . (15) 

From Fig. 1 we can see that the DSNB fluxes being redshifted, arrive on Earth predominantly 
within the energy window Ey = (0 — 35) McV. above which the fluxes are negligible. In this energy 
range water Cerenkov detectors also register events coming from a myriad of other sources. The 
main sources of particles which would imitate the DSNB signal include reactor Ug, atmospheric 
and Pg, solar i/g, spallation products induced by cosmic ray muons, and neutrinos from "invisible 
muons" produced by atmospheric and P^. These form a background for the DSNB signal. 
Events due to reactor z/g appear roughly in the energy range (1.8 — 8) McV and these events can 
be estimated using the information from reactor power and their distances from the detector. In 
the case of SK for instance, it will be even easier to estimate them since KamLAND [104] directly 
observes these events. The events due to atmospheric Ug and Pg are expected to be lower compared 
to those due to DSNB below £^ ~ 30 McV. Number of events expected from atmospheric z/g and 
Pg can be anyway estimated using the predicted fluxes at these energies and can be included in the 
analysis of DSNB events. Therefore, these events do not pose a very serious threat to the DSNB 
analysis. Events due to neutrinos coming from the Sun fall in the energy range E^, ^ 20 MeV and 
can also be estimated fairly well using the fluxes from the standard solar model as well as from 
the direct observation of the *^B fluxes at SNO [105] . These neutrinos can also be identifled in the 
detector from their directionality. Indeed these arc the solar neutrino events that SK observes. 
Therefore, these events do not pose a serious threat to DSNB observation either. The type of 
events which cause a serious concern are the ones produced from spallation. These events are 
typically important in the energy window relevant for solar neutrinos, viz. for ^ 20 MeV. The 
SK collaboration in their paper [67] show that after suitable cuts there are almost no spallation 
events above Eg > 18 MeV. The lower threshold for the neutrino energy is hence restricted to 
El, > 19.3 MeV. The upper limit is taken as 30 MeV. 

Despite the different cuts and selection criteria there are two sources of neutrinos which still 
appear as backgrounds for the DSNB detection. The first has already been discussed above - the 
Ug and Ug events from atmospheric neutrinos. These background events have to be estimated using 
the detector Monte Carlo. The second type of background comes from "invisible muons" produced 
by atmospheric z/^ and u^. These are events where atmospheric z/^ and/or z/^ produce muons with 
kinetic energy less than 53 MeV, which is the threshold for emitting Cerenkov photons. These 
muons/antimuons therefore pass undetected and eventually decay into electrons/positrons which 
are observed by the detector. Estimates for the background due to both these sources have been 
made by the SK collaboration and can be found in [67] . 

2.3.2 Liquid Scintillator Detectors 

Number of events expected in liquid scintillator are also given by Eq. (13). The predominant 
reaction is Ug capture of protons (cf. reaction (12)). The other detection reactions in liquid 
scintillators are charged and neutral current scattering off electrons, charged current capture of Ug 
and Pg on ^^C, neutral current break-up of ^^C (see [106] for reactions of ^^C)^ and neutral current 

^Liquid scintillator detectors can also detect the DSNB i^e flux by their charged current interactions on ^^C [62]. 
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scattering off protons [107]. However, the cross-section for these processes is small, especially at 
low energies [62, 108, 109], and we reiterate that due to redshift the DSNB fluxes are peaked at 
lower energies. Therefore, even for liquid scintillators the main detection weapon is the reaction 
(12). However, compared to the water detectors, liquid scintillators can use the reaction (12) more 
efficiently, whereby they tag the released neutron. While the positron is detected promptly, the 
neutron is captured by a proton in the detector, releasing a 2.2 MeV photon which is detected 
in delayed coincidence after 180 //s. This results in lesser problems with backgrounds, and hquid 
scintillator detectors can be used to observe the DSNB neutrino in the broader energy window of 
= (10-25) MeV [55]. 

The other major difference between the liquid scintillator detector and water Cerenkov detector 
is in the energy resolution, which is much better for the former. The HWHM for liquid scintillator 
detectors is expected to be better than 

as = 0.1^£;^(MeV) - 0.8 . (16) 

The KamLAND detector in Japan [104] and Borexino in Italy [110] are the currently running 
liquid scintillator detectors. While KamLAND has a total mass of 1 kton, Borexino is much 
smaller and comprises of about 300 ton of liquid scintillator. The detectors for the upcoming 
second generation reactor experiments designed to probe ^^13 would be far too small to contribute 
to the study of DSNB neutrinos. However, one could look forward to proposals such as LENA 
[55] which would be situated in the Pyhasalmi mine in Finland and is expected to have 50 kton 
of hquid scintillator. Such a big hquid scintillator detector could collect sizable number of DSNB 
events and prove to be a pivotal player in this game. Another large liquid scintillator detector 
proposal is the Hanohano project in Hawaii [111]. 

2.3.3 Gadolinium Loaded Water Cerenkov Detectors 

The neutron released in the reaction (12) when captured on protons emits only a 2.2 MeV photon. 
This is below the detection threshold of water detectors and hence they cannot normally tag the 
released neutron by delayed coincidence, as liquid scintillators can. However, things could change 
dramatically if GdCls is dissolved into the water. Gadolinium has a large cross-section for neutron 
capture and the capture of neutron on Gadolinium releases a 8 McV gamma cascade. This being 
above the energy threshold, could be easy to observe in water detectors [69] , transforming them into 
giant Ue detectors with statistics many times the statistics expected in scintillator detectors. This 
could give exceptional sensitivity to neutrino oscillation parameters using reactor antineutrinos 
[112]. This will help also in DSNB detection by lowering the lower energy threshold, and we should 
be able to use the same energy window as in liquid scintillators. Following [55], we present our 
results for the energy range (10 — 30) MeV. The energy resolution of course continues to be given 
byEq. (15). 

2.3.4 Liquid Argon Detectors 

Liquid argon TPCs are unique as they allow the detection of Ug. The only other sensitive 
detector technology that we have so far seen built on a large scale was the heavy water detector 
at SNO. However, SNO is now dismantled. Significant amount of R&D on the other hand has 
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gone into the liquid argon option. The ICARUS detector [113] in Italy already consists of a 600 
ton module and has shown the feasibility of this detector technology. Since it is one of the few- 
detector types which can be built on a large scale and allows for very fine granularity, good electron 
detection efficiency as well as detection of r events, this is often considered as a far detector option 
for the Neutrino Factory. Feasibility of probing galactic SN neutrinos was studied in [114, 115] and 
DSNB in [63, 62]. A future large liquid argon detector could have a mass of about 100 kton. Some 
of the currently pursued proposals include GLACIER [116], MODULAr [118] and FLARE [117]. 
Energy resolution in this detector is expected to be extremely good and at energies relevant for the 
neutrino factory it is believed to be in the ballpark of (Te ~ 0.03^£' (GeV) Therefore, at energies 
relevant for DSNB, one can assume that the energy reconstruction could be almost perfect. In 
what follows, we work under this assumption and give our results in terms of the neutrino energy. 
Since further R&D would be needed to determine the backgrounds in this detector, we will show 
results for an energy window of (20 — 40) MeV. 

3 Expected Events from DSNB 

3.1 DSNB Antineutrino Events in Water and Scintillator Detectors 

We give in this subsection the number of DSNB Pg events expected in water and scintillator 
detectors. The total expected number of events are presented in Table 3. The energy windows in 
which we have calculated the total number of events were discussed in the previous section and 
are shown in the parentheses in the first row of the tabic. The number of events per year have 
been calculated assuming a fiducial mass of 22.5 kton for SK and Gadolinium loaded SK (GDSK), 
1 Mton for a future megaton detector (marked in the table symbolically as HK) and Gadolnium 
loaded megaton water detector (GDHK) and 50 kton for the scintillator detector (LENA). The 
results for NH remain the same for any value of ^13. For IH the neutrino oscillation probability 
and hence the number of events depend on 9i^,. We explicitly show results for two extreme values 
of - for small 61^ such that the jump probability P13 = 1 and for large 6*13 such that the 
jump probability P13 = 0. For showcasing the impact of collective effects on the predictions for 
DSNB (anti) neutrino events, we also present in the Table 3 expected number of events if collective 
effects were not taken into account. These arc shown in parenthesis. When there are no collective 
effects, one has only standard MSW transitions in the SN and it is well known that in this case 
antineutrinos undergo maximal flavor transitions for IH when ^13 is large (P13 = 0), while for small 
values of ^13 or with NH (for any ^13) there is no matter enhanced resonant oscillations and these 
two scenarios give identical results. We therefore get larger number of events for IH and large ^13. 
However, once the collective effects are switched on, the small and large ^13 cases of IH switch 
roles. Since there are now two stages of ffavor conversions, first due to collective effects deep inside 
the SN and then due to MSW transitions, the final fiuxes are such that IH with large 6*13 and 
NH give identical predictions, while IH with small ^13 predicts larger number of events (cf. upper 
right panel of Fig. 1). 

It can be seen, that we expect about a couple of events per year in SK ^. This would go up by a 

^Note that there will also be a large number of background events in the detector and one has to find the signal 
by looking at excess of events above the fluctuations in the background. This makes DSNB detection more difficult. 
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6.2 (6.2) 


98.2 (98.2) 
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7.7 (7.7) 


LL 


IH (Pi3 = 0) 


2.5 (4.4) 
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98.2 (175.7) 
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IH (Pi3 = 1) 


4.4 (2.5) 
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Table 3: Number of expected events per year per 22.5 kton of SK and GDSK, 1000 kton of HK 
and GDHK, and 50 kton of LENA. The events without collective effects are shown in parenthesis 
for comparison. 
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Figure 2: Number of expected events as a function of the jump probability Prs for the Gl model. 
Black lines are for NH and blue dashed lines for IH. The yellow dashed dotted lines show the case 
for IH without collective effects (WOC). 

factor of about 2 — 3 if Gadolinium were to be added to the water. The corresponding number for a 
megaton of water would be scaled upwards by a factor of 1000/22.5 and we expect about 40 — 176 
(140 — 356) events per year in megaton water (Gadolinium loaded water) detectors depending on 
the choice of the neutrino mass hierarchy and 6'i3 and the SN model. After 10 years of running 
these numbers would be a factor of 10 higher, and we could have a few thousand events in the 
Gadolinium loaded detector. It should therefore be straightforward for megaton water detectors, 
with or without Gadolinium, to be able to observe these DSNB fluxes. More importantly we note 
that for a given SN flux model, it should be easy for megaton water detectors to determine the 
hierarchy, if sin^ 6'i3 ^ 10^^. For almost vanishing ^13, we can see that for Gl, NH predicts 1960±44 
(678 ± 26) events in 10 years of running of GDHK (HK) while IH predicts 2960 ± 54 (1096 ± 33). 
It would therefore be easy to distinguish one hierarchy from the other. Note that this is one of 
the very rare type of experiments which can give information about the neutrino mass hierarchy 
even if ^13 was below the reach of the most Neutrino Factory and Beta-beam experiments. A 50 
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Figure 3: Number of expected events per year in 2 MeV positron energy bins in SK (upper panel) 
HK (middle panel) and LENA (lower panel). The solid black lines show the projected event 
spectrum for NH while the dashed blue lines are if IH was true with P13 = 1. The SN flux model 
corresponds to Gl. 

kton liquid scintillator detector should be able to record 46 — 106 events in 10 years of running. 

We have shown in the Table, number of events expected assuming either the Gl, G2 or LL 
model for the initial SN neutrino fluxes. We find that the lowest number of events are predicted 
by the G2 model, while LL predicts the highest event rate. In fact, one can see that the event rate 
predicted by NH and Gl is close is that predicted by IH and G2. Likewise, the rate predicted by 
IH and Gl is close to the one predicted by NH and LL. We have discussed before the uncertainty 
associated with the SN models. Therefore, if the uncertainty in the model predictions for the 
initial fluxes remain at the current level, then it might be hard to distinguish the hierarchy from 
the DSNB itself, especially in the smaller detectors. However, for Gadolinium loaded megaton 
water detectors it might still be possible to say something about the hierarchy. Also, for G2 and 
NH (LL and IH) we have a prediction which is lower (higher) than any other case and therefore for 
these cases there is no confusion. For instance, if GDHK records less (greater) than 1500 (3000) 
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Gl 


G2 


LL 


NH fPiQ = 0) 


4.9 (4.9) 


2.3 (2.3) 


9.9 (9.9) 


NH (Pi3 = 1) 


3.6 (3.6) 


1.7 (1.7) 


7.3 (7.3) 


IH 


4.9 (3.6) 


2.3 (1.7) 


9.9 (7.3) 



Table 4: Number of Vf. charged current events on ^°Ar per year per 100 kton of Liquid argon TPC 
in the energy window — (20 — 40) MeV. The events without collective effects are shown in 
parenthesis for comparison. 

events, we could say that the hierarchy is normal (inverted). Of course, we have nowhere taken 
into account the uncertainty in the star formation rate. That might bring additional complication, 
which we do not address in this paper. 

So far we have presented results only for two extreme cases of ^13, very low corresponding to 
Pi3 = 1 and very high corresponding to P13 = 0. For intermediate values of the mixing angle 
the jump probabihty ranges between and 1. We show in Fig. 2 how the total event rate in the 
different detectors change as a function of P13. The SN model assumed is Gl. Solid black lines 
show the case for NH while the dashed blue lines show the case for IH, where we have included 
both collective as well as MSW transitions inside the SN. It is easy to see from the expressions 
given in Table 2 that for IH, the event rate would rise almost hnearly with P13. If collective effects 
were not taken into account then the trend would have been the opposite, and we would see a 
decrease in the v^. event rate with P13. These are shown for the different detectors by the yellow 
dot-dashed lines in the figure. 

For sizable number of events, it might even be possible to do a spectral analysis of the DSNB 
events. We show in Fig. 3 the event spectrum for 22.5 kton SK (upper panel), 1 Mton HK (middle) 
and 50 kton LENA (lower panel). The events per year are shown in 2 MeV energy bins. The solid 
black lines give the event spectrum for NH while the blue dashed lines are for IH with P13 = 1. We 
show results where both collective as well as MSW oscillations are taken into account. Upper and 
lower energy threshold for the different cases are indicated by vertical lines and we have assumed 
the Gl model for the initial fluxes. 

3.2 DSNB Neutrino Events in Liquid Argon Detectors 

Liquid argon TPC could offer a unique laboratory to probe from a future galactic SN as well as 
from the DSNB around us. We show in Table 4 the number of expected v^. charged current events 
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Figure 4: Number of expected events per year in 2 MeV neutrino energy bins in a 100 kton Liquid 
argon TPC The solid black lines show the projected event spectrum for NH (P13 = 1) while the 
dashed blue lines are if IH was true. The SN flux model corresponds to Gl. 



on ^°Ar. We show results for NH and large (P13 = 0), NH and small ^13 (P13 = 1), and IH 
for any value of ^^13. Expected number of events are shown for the three benchmark flux models. 
We see that the number of z/g events expected in liquid argon detectors is extremely small. This 
is because the z/g +'^° Ar — »• e~ K* cross-section (taken from [115]) is very small at low energies 
and rises very fast as the energy increases. The DSNB flux on the other hand gets redshifted to 
lower energies thereby reducing the number of events. In particular, the DSNB flux is peaked at 
around 5 MeV, with very few neutrinos in the energy window 20-40 MeV (cf. Fig. 1). It might 
also be interesting to compare the number of expected z/g events in a liquid argon detector, with 
the number of Ve events in a water detector, for same number of target nuclei/nucleons. It turns 
out that 100 kton of liquid argon has 1.5 x 10^^ argon targets, while 22.5 kton water detector (SK) 
also has 1.5 x 10^^ proton targets. On the other hand, the cumulative cross-section in the energy 
window of 20-40 MeV for z/g capture on ^''Ar is larger than the cross-section for z/g capture on 
protons by a factor of about 2. Signal in this energy window, for the LL SN model with complete 
flavor conversion^, would be 9.9 and 5.8 events in 100 kton of liquid argon and 22.5 kton of water, 
respectively. This implies a ratio of about 1.7, which agrees with the rough estimate of the factor 
of 2 coming from the difference in the cross-sections. If we could lower the energy threshold in 
liquid argon to 5.5 MeV, we could expect about 8.1 events per year for NH with small di^ and 
about 10.5 events per year for IH (and NH with large ^13). In Fig. 4 we show the event spectrum 
in bins of 2 MeV width. The black solid line shows the spectrum for NH with P13 = 1 while the 
blue dashed line is for IH. 

^For complete flavor conversion, the resultant flux at both liquid argon and water detector is Vx, and is therefore 
the same. Of course in reality, complete flavor conversion can be possible only in one channel. The above example 
is just to illustrate the difference in the number of events for the two detector types being compared. 
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4 Summary and Conclusions 



Neutrinos emitted by core-collapse supernovae over the entire history of the Universe, pervade us. 
This is the so-called diffuse supernova neutrino background. The DSNB fluxes are theoretically 
given by folding the neutrinos emitted from a typical SN with the rate of SN explosions as a 
function of the redshift, and integrating over all redshifts to take into account all possible SN 
explosions that might have happened in the Universe. Since the final fluxes emerging from the SN 
depend on neutrino flavor coversions inside the SN, the DSNB fluxes also depend very crucially on 
neutrino properties. Therefore, while a galactic SN event is eagerly awaited in order to shed light 
on SN theory on one hand and neutrino properties on the other, detecting the DSNB in currently 
running and future detectors could help us constrain SN dynamics, cosmic star formation rate 
as well as neutrino properties. However, the primary agenda is to succesfuUy observe them in 
terrestrial detectors. 

Being redshifted, the spectrum of DSNB fluxes is peaked at smaller energies, making their 
detection even more challenging. So far the running Super-Kamiokande detector has managed to 
put an upper bound on the Pg DSNB flux. However, the situation might improve in the future with 
possibility of a signal in the upcoming large scale detectors which would be built to serve as the 
far detector for high performance neutrino beam experiments. Observing DSNB would be free for 
these detectors and the physics output from that would be immense. In this paper we re-analyzed 
the potential of a selected class of future detectors to detect DSNB fluxes. Such an analysis has 
been warranted by the flurry of activity in the held of SN neutrino research, following the revival 
of interest in neutrino-neutrino self-interaction. These interactions inside the SN have been shown 
to produce significant change to the final neutrino fiux spectrum, especially if the neutrino mass 
hierarchy is inverted. Since these so-called collective effects inside the SN are unavoidable, it was 
necessary to revisit the issue of DSNB detection. 

We considered water. Gadolinium loaded water and liquid scintillator detectors for DSNB 
detection and liquid argon TPC for observing the DSNB flux. A major issue in this held is 
the model uncertainties in the SN neutrino fluxes themselves. We presented results for three SN 
neutrino flux models. We calculated the total number of events for both the hierarchies and for 
two extreme values of ^^13 resulting in jump probability P13 and 1. Number of events expected 
in future megaton water and 50 kton liquid scintillator detectors are large, with a few thousand 
events expected in Gadolinium loaded megaton water detectors running for 10 years. For true 
inverted hierarchy, it becomes possible to get very large flavor oscillations even if ^13 — > 0. We 
showed that under fortunate circumstances, it might be possible to get information on the neutrino 
mass hierarchy by observing DSNB in megaton water detectors. Note that this is a very unique 
situation, since for 6'i3 — > it becomes almost impossible to determine the hierarchy using long 
baseline experiments. In this way, DSNB detection could be complementary to the long baseline 
program. We also showed how the total number of events change if ^13 increases from very small 
to very large values, decreasing P13 from 1 to 0. Finally, we showed the event spectrum by binning 
the prospective data in 2 MeV bins. 

In conclusion, very large number of DSNB events are expected in the next generation detectors 
and therefore, it should be possible to observe DSNB in the future. Collective effects inside 
SN signiflcantly change the predicted number of DSNB events if the hierarchy is inverted. Under 
fortunate conditions it might be possible to determine the neutrino mass hierarchy using ths DSNB 
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signal and this could be done even if 6*13 0, in which case long baseline experiments would not 
be able to tell the hierarchy at all. Propects of DSNB detection look extremely promising and 
one might even feel optimistic about learning about neutrino oscillation parameters, cosmic star 
formation rate and maybe about SN physics, by observing these relic neutrinos in future detectors. 
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